On resolution to Wu's Conjecture 

by Theodore Yaotsu Wiu 

Abstract. In this series of studies on Cauchy's function f(z) {z = x + iy) and its integral J[f{z)] = 
(27ri)~^ <f^ f{t)dt/{t — z) taken along a Jordan contour C, the aim is to investigate their comprehensive 
properties over the entire z-plane consisted of the simply-connected closed domain T)^ bounded by C 
and the open domain T)~ outside C. This article attempts to solve an inverse problem that Cauchy 
function f{z), regular in T)'^ and on C, has a singularity distribution in 'D~ which can be determined 
in analytical form in terms of the values f{t) numerically prescribed on C, which is Wu's conjecture[l]. 
^) \ It is resolved here for f{z) having (i) a single, (ii) double, or (iii) multiple singularities of the types 
C5 ■ (I) T,j^Mj{zj — z)^J , (II) M£log{z2 — z), by having their power series expanded in z and matched 
on a unit circle {t = e*^, — vr < 9 < ir for contour C) with the numerically prescribed Fourier series 
f{z) = Sg°c„e*"^ for solution. The mathematical methods used include (a) complex algebra for cases 
I (i)-(ii), (b) for case (iii) a general asymptotic method developed here for resolution to the Conjecture by 
' induction, and (c) the generalized Hilbert transforms to expound essential singularities. This Conjecture 
has an advanced version for f{z) to be given only one of its two conjugate functions on C to suffice, and 
another for the complement function F{z) defined as being regular in domain P~ and having singularities 
in V~^. These new methods are applicable to all relevant problems in mathematics, engineering and 
mathematical physics requiring breakthrough by having the exterior singularities resolved. 
Key words: Cauchy function, singularity distribution, Wu's Conjecture, Resolution by induction. 
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1. Introduction. 

The foundation of the present study starts with Cauchy's theorem and integral formula. 



eg : Ami - ^ / 



J[f(.z)] = — / -^^dt = f{z) {z£V+ ~ open domain inside C), (la) 

2vri Jq t — z 

^ I = (z € — open domain outside C), (lb) 

I where Cauchy's function f{z) is assumed to be analytic, regular Vz G D'^ and continuous for z = t on 
' contour C, taken in positive (counter-clockwise) sense as understood. Here, (jlbp follows from Cauchy's 
. integral theorem that g{t)dt = if g{t) is regular within and on contour C (as is g{t) = f[t) /{t—z) Vz € 
P~), while (jlap is known as Cauchy's integral formula, also called Cauchy's functional relation, known 
^ to hold for z in open domain T>^ (but not including its boundary contour C). 

^ • The task of determining the value of Cauchy's integral J[f{z)\ for z situated right on C has been 

^ ■ accomplished by Wu[l] with adopting a generalized condition that 

f{z) be C" Vz G and in a neighborhood Mc striding across contour C (n being arbitrary), (2) 

where the corresponding function f[z) is called the generalized Cauchy's function. This new condition 
thus renders the integral J[f{z)\ = {27ri)^^ f{t)dt/{t — z) intact in value as a point z G T)^ (or 
z € T>~) tends to a generic point Zo on C without crossing the contour if the contour remains fixed 
except for an infinitesimal stretch about Zq where it is indented into a semi-circle onto the P~-side 
(or a semi-circle C~ onto the 'D"'"-side), of radius e from Zq on C (see Wu[l], Fig. 1). In the limit as 
e — s- 0, the integral J[f{z)] having zjs P+j Zq (from inside C) is let to assume its inner limit of 
undetermined value f~^{zo) (resulting from the integral over the semi-circle plus the integral over 
C—C^ which assumes its principal value as e — > 0). Similarly, the integral J[f{z)] having z{£ T)~} Zq 
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(from outside C) has its outer limit of J[/(zo)] = f~{z) = as is invoked by Cauchy's Theorem (jlbp . 
and this hmit is the sum of the integral over the semi-circle C~ plus the principal value of the integral 
over C — C~ . The final two limit equations thus yield three key relations (cf. Wu[l], Eq. (12)) as 

(!) : /+(z) = fiz); {II) : /-(z) = 0; (///) : /(z) = -V <f j^dt {z G C), (3) 

m Jet- z 

in which the integral with symbol V denotes its Cauchy principal value while the suffix of zq is omitted 
for all z on C. Here, relation (I), f~^{z) = f{z), shows that the limit f~^{z) of integral J[f{z)] resulting as 
z {€ T>~^} — > Zq is equal to the original prescribed f{z) \/z G C, therefore proves the uniform continuity 
of f{z) in the dosed domain T>+ = [D^ + C], relation (II) shows J[f{z)] = V z in the closed domain 
T>^ = \D~ + C], whilst relation (III) relates f{z) for each z C in terms of all the other values of f{t) 
over C. Since the original Cauchy's two integral theorems in ([T]) are thus shown also holding valid for 
z G C, this further establishes the theorem that Cauchy's integral J[f{z)] is uniformly convergent in 
the closed domain D+ = [P+ + C]. These are among the important results discovered in Wu[l]. 

The key relations (I)-(III) in ([3]) have various prospects for further applications and development. 
When the contour C, originally arbitrary in shape, assumes certain geometric forms in particular, e.g. 
one that circumventing the upper-half (or lower-half) of the z-plane, or another inside (or outside) the 
unit circle \z\ = 1, the specific Hilbert transform relations between the conjugate functions u and v of 
f{z) = u{x,y) + iv{x,y) follow directly from the formulas in ([3|) (cf. Wu[l]), which can therefore be 
regarded as providing the generalized Hilbert transform formulas with a simple proof. 

The relations (I)- (III) further help illustrate the fact that a unique relation exists between the values 
of an analytic function f{t) \/z = t varying on contour C and all the singularities of f{z) \/z G T)~ outside 
C, as exemplified in explicit cases of the so-called direct problems defined by having the singularities 
prescribed (cf. Wu[l], §7). The only exception to this assertion of ever existence of singularities of 
f{z) Mz G T>^ is when f[t) = 1 (or any constant) on C, in which case f{z) = 1 Vz in the entire z-plane, 
including z = oo, as can be shown by analytic continuation, in virtue of Liouville's theorem. 

The present study is motivated by Wu's[l] inverse problem asserting that every Cauchy function 
f{z) invariably has a singularity distribution in domain T>^ outside C which can be determined in closed 
analytical form in terms of the values f{t) prescribed only numerically on contour C, i.e. 
The inverse problem. The inverse problem is (i) to have function f{t) prescribed over contour C 
only numerically, or equivalently in terms of a series with known numerical coefficients, (ii) f{z) is 
regular inside contour C , and (Hi) to use the given numerical data to determine the entire singularity 
distribution of f{z) Vz G T>^ outside C in a closed analytical form, whatever the singularity distribution. 
This seems to cover such unsettled situations as resulting from using, e.g. perturbation series expansion, 
boundary integral scheme, or other computational algorithms. For the general solution to the inverse 
problem, its existence is enunciated in terms of the following Conjecture (cf. Wu[l], §10.3): 
The conjecture. A solution to this inverse problem is conjectured to exist, which may not be unique. 

Resolution to this inverse problem is of vital importance. This can be well expounded by an 
expository review of the interesting and colorful history of the theory of water waves, a subject of 
intimate significance to, and long desired for the general resolution to the Conjecture at hand. 

Historically, it dates back to the pioneering work of Sir George Gabriel Stokes [2] who developed the 
perturbation expansion theory (1847) and used it for the very first time to study the nonlinear effects in 
simple harmonic water waves on deep water to three terms of the power series (expanded in the physical 
plane) in terms of a small parameter a = ka, a being the wave amplitude, k the wave number. The 
crucial results attained there were found to produce important qualitative changes in wave behavior and 
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manifest new phenomena, such as the dispersion relation involves the wave amplitude. Later in 1880, 
Stokes pursued the expansion for periodic waves on water of constant rest depth h, this time formulated 
in the complex potential plane (for its prefered fixed boundary of the flow region, thereby lending much 
simpler algebra without involving iterations) to five higher order terms [2]. Along this line, perturbation 
expansions have been carried out for solitary waves by various authors, including Longuet-Higgins &: 
Fox [3], Williams [4] and others to further higher orders as summed in review [5], yet still falling short of 
a conclusion on the series convergence. Recently, the series expansion for solitary waves of all heights 
have been brought by Qu[6] to fifteenth order in a series of exact coefficients in rational numbers. This 
has been further pursued by Wu et al.[7] to discover that the series is asymptotic in nature within the 
range of 12 — 17th order (varying due to which one or another small parameter adopted for expansion). 
But for all the first integrals (e.g. of the excess mass, wave energies, etc), convergence was found by 
Wang & Wu[8] to arise within the range of 20 — 25th order of the power series expanded only in terms 
of e = k'^h^, a final result which is not shared by power series in other parameters, e.g. a/h, attested. 

On the other hand, the issue on the exterior singularities has also been apprehended as a new 
interest underlying this general problem. Grant [9] considered the highest periodic water wave (with 
Stokes's 120° corner crest), formulated with the physical variable as a function z = x + iy = z{f) 
in the lower half of the complex potential f = (f) + iijj plane continued analytically to the entire /- 
plane. He showed by analysis that the singularities of z{f) outside the flow field 2?+ are all of the 
regular type of order 1/2, i.e. z{f) c(f — foy^"^ sls f ^ fo ^fo ^ 1^^- Grant contended that there 
exist several singularities of order 1/2 ^ 2?+ which coalesce at the cornered wave crest to become 
a singularity of order 2/3 plus a secondary singularity of an irrational order, hence not regular, in a 
quite perplexing way. Along a computational approach, Schwartz [10] examined the singularity using 
the Fade approximation and Domb-Sykes[ll] plots, flnding that the singularity varies continuously with 
the wave height from order 1/2 to 1/3 (for the complex velocity df /dz), and drawing a conjecture that 
the changes involve coalescence of several singularities of order 1/2. However, Tanveer[12] examined 
the number of singularities, based on the argument principle, finding that there is just one singularity 
outside the flow field, in disagreement with Schwartz's results. This question was further investigated by 
Longuet-Higgins &: Fox[13] with analytic continuation of their z{f) for almost-highest wave expanded 
in a power series to sixty terms with the Fade approximation to find an exterior branch cut of order 
1/2, in agreement with Grant[9]. However, the conjectured coalescence of singularities into one to fit 
the cornered crest of the highest wave remains to be expounded. Finally, any exterior singularity must 
have its location and magnitude, in addition to its type, also determined to complete the resolution. 
This is the main objective of the present study. 

Summing up the foregoing two issues, it appears that their eventual resolution may be attained 
together with first resolving the Conjecture at hand. 

We now return to our present endeavor for resolution of Wu's Conjecture. First, we classify function 
f{z) to have (i) a single, (ii) double, or (iii) multiple singularities of the algebraic and/or logarithmic 
types, each of magnitude Mj, located at Zj, of type kj, to be determined in terms of the numerically 
prescribed f{z) at the unit circle z = e*^ as the standard contour C. For cases (i)-(ii), solutions can 
be obtained using complex algebra for each singularity of f{z) in outside C, as illustrated in §2- 
§3. The original Conjecture is extended to an advanced version for f{z) to have only one of its two 
conjugate functions prescribed on C to suffice, as shown in §5, and further for the complement function 
F(z) defined as being regular in domain 'D~ and having singularities in P"*". The case for singularities 
located on a circle in 2?" concentric to C is discussed in §6. For case (iii), solutions by using complex 
analysis are shown in §7, and also by a general method which is developed in §8 for resolution to the 
Conjecture by induction. Finally, the essential singularities are expounded in §9. 
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2. Generalized Cauchy's function with a single singularity. 

First, we consider the primary case when f{z) (regular within and on contour C) possesses a single 
singularity (of number Ng = 1) to be determined in terms of its values f{t) prescribed numerically for 
z = tMt ^ C , however with no further data available on the base parameters (M, zi , k) for magnitude 
M, location zi, and primary key index k (for the order, e.g. M(zi — z)^) of the sole singularity. While 
being primary, this is the basic case to which more general cases can be reduced or referenced. 

For simplicity, we take the unit circle |t| = 1 for the standard contour C throughout, onto which 
any arbitrary regular contour enclosing a simply-connected domain can be mapped by conformal trans- 
formation, so that in polar form, z = rexp{i9), r = \z\, = argz, r = 1 on C, on which we have 

f{t) = f{^') = f{e) i-7r<e<7T), (4) 

which is prescribed numerically to a desired accuracy so as to satisfy condition ([2]). 

With this provision, we can then express f{6) in a Fourier series in complex form as 



N 

J{e) = cne*"', ^ = ^ f e-'^'f{e)de, (5) 



n=0 



with N numerically set for the accuracy desired. This series representation may be regarded as entirely 
equivalent to f{9) given numerically in (jH), for they are merely a quadrature apart. 

On the other hand, the sole point singularity can be represented generically as 

f{z) = M{zi-zf (ki|>l), (6a) 

h{z) = M,log(zi-z), (6b) 

where k R, assumed real, so that it is a zero (or a pole) for k being a positive (or a negative) integer, 
or an algebraic branch connecting zi and z = oo for k not an integer, or a logarithmic branch for A; = 0, 
all lying outside the unit circle for \zi\ > 1. Thus, (|6a p - ()6bp represent f{z) with a single singularity for 
resolution, i.e. to have {M,zi,k) determined in terms of the given coefficients {c„}. 
To proceed, we also expand f{z) of ^ into a power series in z, yielding 

^ fcV- f fk\ k{k-l)---{k-n + l) r{k + l) 

1 / z \ " 

/^(z) =M^log(zi-z) = M^log(2i)-M,V- - , (7b) 

where r(A;) is the Gamma function, r(A; -|- 1) = kT{k), T{n -|- 1) = n!. 
Next, matching series (I7ap termwise with series ([5]) readily yields 



/(z) = L^7„, ( -— j =s{z), Co = L70, Cn = Ljni-Zl) (L 



Mzf) 





agreeing with (0) for z = e*^ on contour C. This series has the ratio of consecutive coefficients as 

Rn{k,zi) = = = — — -— = 1 —r — (n = l,2,---), (8b) 

Cn In zi n + lzi \ n + lj Zi 
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which shows (e.g. by the ratio test) that series ([7|) converges for \z/zi\ < 1 within the circle of con- 
vergence of radius p = \zi\ on which the point singularity lies at z = zi. Finally, the undetermined 
singularity location zi can be eliminated by taking the ratio Dn(k) = Rn+i/Rn, giving 



Rn+1 n + ln + 1 — k 



Dn{k) = = ^ _ ^ = Cn+2 c„/c„+i (n = 1, 2, • • • ), (8c) 

from which k is attained in terms of the given data Cn+2Cn/ c^+i to identify the singularity type. 

Further, it is clear that the preceding relations Q-dS]) for f{z) also hold valid for the sole logarithmic 
singularity f^{z) of pbj) for A; = 0. 

At last, we note that if (i) f{z) has a single singularity, (ii) the given numerical data {c„}'s are 
exact, then the above relations for n > 3 all turn out to be identical in effect to the first three orders 
for n = 1,2,3 which can provide the final solution. We therefore have the following theorem proved. 

Theorem 1. Primary resolution to the Conjecture. If (i) function f{z) is regular in domain T>~^ 
bounded by the unit circle C {z = e*^) and on C, (ii) f{z) is numerically prescribed on C by series ^ 
with coefficients {cn, n = 0, 1, 2, • • • } which are exact in value, and (Hi) f{z) has a single singularity of 
unknown magnitude M , location z\ i\z\ \ > 1), and key index k (of order {z — zi)^), then 

' = "+ n+l-(t + 2)Z)„ ''" = ^' ''" = ^' 

^1 = t l^fg (=^iN)' M = co/z'l, M^ = co/log(zi). (9) 

In principle, the formulas in ([6|)-([9|), if exact, should produce equal values to {M,zi,k) for every 
order n. In practice, however, the numerical data for the coefficients {c„}'s may contain errors. Even 
more so, the presumption for the present case of A'^s = 1 that there is only one singularity needs first to 
be confirmed. In this premise, it raises the following important issue to be addressed next. 
3. Verification of consistency and error estimate. 

If the final results for base parameters (M, zi, k) obtained from a specific trial run turn out to be 
uniformly consistent over order n except for very small fluctuations from order to order, which may be 
due to numerical errors incurred in deriving the coefficients {c^j's for prescribing the given data. In 
such cases, the resolution may be said to bear the confirming feature of consistency in verifying the sole 
hidden singularity. If in contrast, the resulting values for (M, zi, k) appear to vary in finite magnitude 
for differing orders, such inconsistency then indicates that the correct solution is yet to be found. 

For error estimate, we assign a small error e„ to each Cn to give the rounded coefficient Cn and 

Rn = Cn+l/Cn = iCn+l+en+l)/{Cn + en), > k{n) = k + K{en, Cn+l, €n+2) , (10) 

SO that /t, the error to k, being by ([Hell an algebraic function of three specific e's, vanishes with the e's 
instead of being nonlinearly magnified. A similar conclusion can be drawn for {zi, M, M^). 

Example 1. To illustrate these points, we let series ([5]) be prescribed for a specific case by 

Co = 2, ci = 1.999895, C2 = 1.333415, C3 = 0.889123, C4 = 0.592593, C5 = 0.395062, 
for the five leading terms taken here for solution and verification. We then find, by ([9|), 

= C2/C1 = 0.666743, R2 = 0.666801, R3 = 0.666492, R4 = 0.666669, 
Di = R2/R1 = 0.999912, D2 = 0.999535, Ds = 1.00026, 
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and by Q, k{n) can be determined for n = 1, 2, 3 as 

k{l) = -0.999823, k{2) = -1.002791, A;(3) = -0.996851. 

As a result, it seems pertinent to regard each of these /c(n)'s retaining such a uniformly small error from 
the round figure of A; = —1 to lend it as the rounded value for k with sufficient verifications. Similarly, 
we may let each of zi{n) = 1/Rn assume the rounded value of zi = 1.5 with supporting checks. We 
therefore obtain the final result as 

k = —1, zi = 1.5 M = ziCq = 3, 

in resolving the Conjecture for this case with f{z) involving just one simple pole at zi = 1.5, of 
magnitude M = 3. Finally, we can draw from the rounded coefficients c„'s their specific errors as 

c„ = c„ + e„ — > ei = 1.05 X 10~^ e2 = -8.2 x 10-^ ea = 2.35 x 10"^, • • • , 

with similar error estimates for [Z)„ — Dn\-, 9 Dn — Dn — > V e„ ^0. Here we point out that the final 
solution can be used to re-generate the Fourier coefficients for an ultimate verification of the errors. 

We also see that this case can suffice to exemplify all the others in this group. In general, it needs 
merely several leading terms of series ([5]) to complete the relevant resolution as so illustrated. 

4. Advanced Conjecture for the generalized function f{z). 

We can extend the original Conjecture to another level of higher interest as follows. 
Advanced Conjecture: The new inverse problem is to have only one of the conjugate functions 
^/ fi^) — ''J'{x,y) + v{x,y), say v numerically prescribed on contour C of the unit circle (on which 
z = e*^,w = v{0) ,—iT < 6 < tt), other premises being equal, then the new Conjecture asserts that the 
exact singularity distribution of f(z)yz G outside C exists and can be determined. 

For resolution, we first separate the series variables into real and imaginary parts as 
Cn = CLn + ibm CnC^"'^ = {on COS nO — bn sm n9) + i{bn COS n9 + Gn s'm uO) = An + iBn, 

N 

f{e) = Cne'"' = Y, ^n + iBn = u{e) + iv{e)- (11a) 

n=0 n=0 

With v{9) given, u{6) can be derived by applying the circular Hilbert transform (cf Wu[l] Eq(33)): 

u{e) = H[vm = ^f cot vi<P) = H-'HO)] = f u{e) cot ^de, (lib) 

which holds for normalized u and v such that u{(j))d(j) = J^^ v{(f))d(l) = 0. Hence 

u{9) = H[v{(f))] = H[smn(l)] = cosnO, v{(f>) = H~^[cosn9] = —smn(f) (n > 1). (He) 

Using these relations with the leading term for n = already known in (jllap yields 

u{e) = H[v{^)] — > f{9) = u{9) + iv{9) {-TT <9< tt), (12) 

thus proving that given v{9) renders f{9) = u{9) + iv{9) completely prescribed over contour C, and 
hence validating the final solution in ([9]) for the primary resolution to the Advanced Conjecture. 
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5. The Complement Conjecture for the complement function F{z). 

The analytical scheme devised above for resolving the Advanced Conjecture can be further extended 
to exhibit the symmetry in having the roles of domain P"*" and 'D~ interchanged by converting the 
Conjecture for the function f{z) reflected into the contour C of unit circle to the Complement Conjecture 
for the complement function F{z), defined as being regular of order C^^z € and in a neighborhood 
striding across contour C (cf. Wu[l]). Being so defined, it starts with its integral theorem that 

J~[F{z)] = j F{z)dz = - j F{z)dz = 0, (13) 

where contour C~ is on the unit circle \z\ = 1, taken in clockwise direction in its own positive sense. 
By analogy, the complement function must have a singularity distribution inside contour C (unless 
F{z) =const. on C), for which fact we have the analogous conjecture as follows. 

The Complement Conjecture: // the complement function F{z) has on contour C i\/t = e*^) one of 
its conjugate functions of F{e^^) = U{9) + iV{6), say V numerically prescribed, which is equal in value 
to v{9) of f{e^^) = u{9) + iv{9) given on C, then the new Conjecture asserts that the exact singularity 
distribution of F{z)\/z E T>^ inside C exists and can be determined. 

For resolution, we note that invoking V{9) = v{9) dictates that on contour C (it = e*^). 



F{e'^) = F{9) = U{9) + iV{9) = -/(e^^) = -u{9) + iv{9), (14a) 
— > U{e) = -u{9), V{9) = v{9) {-TT<9< tt). (14b) 



where f{t) denotes the complex conjugate of f{t). This therefore implies that on C (Vt = e ), 

- F{t) = m= lie-'') = IS„=o(-^i)-"7n(A:)e-*"^ (= S„=oc„e-^"^) {-tt < 9 < tt), (15a) 

in which the first step follows from ([7]) whereas the second from ([5]), with L = Mz\. This series is 
convergent and can be continued analytically to all z = re*^ over the entire z-plane, thus giving 

- F{z) = Sc„(re*^)-" = Sc„z-" = LS„=o(-^i)""7n(fc)^~" = M (zi - z-')\ (15b) 

where F{z) has its imaginary part equal to that of f{z) on C. This F{z) has a singularity at z = l/zi 
plus a coupled singularity at z = (in refiection to that at z = oo for f{z)) which can be made 
single- valued when k is not an integer with a branch cut between z = and z = 1/zi lying within C. 

Concluding on these extended cases, we note that in addition to the feature that on contour \z\ = 1, 
f{z) and F{z) have their imaginary parts equal and their real parts only opposite in sign, they further 
satisfy Cauchy's two integral theorems as expected, i.e. 

J{f{z)] = — (f ^dt = f{z) {zeV+ = V+ + C), (16a) 



2TTi 



c 



l-Ki 



c 









= V+ + C), 









= v- + c), 







G V+ 


= V+ + C), 


F{z) 






= V~ + C). 



(16b) 



J-[Fiz)] = -^(f^dt = {z eV+ = V+ + C), (16c) 



(16d) 



These integral formulas are presented specifically to hold for each closed domain in virtue of their 
corresponding theorems presented in Wu[l]. 
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6. Resolution for function with complex conjugate pairs of singularities. 

A case of special interest is for function f{z) having a complex conjugate pair of singularities. In 
this case, the series for f{z), in virtue of Schwarz's symmetry, has its coefficients all real. Here it is apt 
to first consider a f{z) having a conjugate pair of logarithmic singularities at zi = e*" and zi, 

f(z) = log[(l - e-'°z)(l - e^"z)]-^/2 ^ 5.00 ^^^n ^ ^.^N a„ = - cos na, (17) 

n 

s{z) being the series function defined by the series (V|z| < 1) versus the sum function f{z) for all z. In 
case when the series function s{z) is the only data available for locating the singularities of f{z) while 
f{z) is still undetermined, this particular series can actually provide a useful basic reference. 
Sign pattern criteria of series, by which certain definite relations can be attained between a (hence 
the argument of zi) and the resulting sign pattern displayed in the series. For this particular series (jl7p . 
the sign of On will depend on the sign of cos(na) and will result in an interlacing sign pattern with first 
A^i terms of positive a„'s, then terms negative, and again terms positive, • • • , so that 

VT TT TT StT 3 1 

--<A^ia<-, -<{Ni + N2)a<—, (j - -)^ < (iVi + • • • + A^i)a < (j - -)7r, 

for j = 1, 2, • • • . In the limit, we actually have just shown a theorem by J. Li[14](1982) in extending the 
pioneering work on this subject by M.D. Van Dyke[15](1974). 

Theorem 2. Li's Theorem on Location of a singularity pair of a series. If the power series 
TiUnz"', with On G R, lim |a„|"'^/"' = 1, has in turn Ni terms positive, N2 terms negative, terms 
positive, ■ ■ ■ , then the series has a conjugate pair of singularities at zi = and zi = e~'", where 

a = lim J7r/(S{^, iV^). (18) 

In particular, ()17p reduces for q = to f{z) = — log(l — z) = Siz", for which A'^i = 00, hence (I18p 
also yields q = 0. Likewise, for a = vr, (fT7|) becomes f{z) = — log(l + z) = Si(— z)", for which 
A^i = = ■ ■ ■ = 1) then j7r/(S-^^-^ Nk) = j'n'/j = vr for j > 2, or as j 00, thus showing the 
singularity located at z = —1 just as that of f{z) = — log(l + z). From these two better known cases, 
we can also infer that a = 7r/2, or a = vr/S, or a = tt/N if the signs of a^'s interlace in pairs, or in 
triplets, or in A^-tuples, respectively. 

In general, the angular positions of all logarithmic singularity pairs can be determined as a sum 
of the individual pairs. Further generalizing Li's Theorem, we can readily resolve the singularities of 
products of multiple pairs of algebraic conjugate functions such as f{z) = {z — zi)~^{z — zi)~'' {zi = e*°), 
since we can then set 

giz) = log f{z) = -A;log{(l - ziz)il-ziz)} = 2k n"^z"cosna, 

which is again in the form of ()17p . and therefore giving the singularities of f{z) as that of f{z) = 
exp{g{z)}. The same logical argument holds for other types of such functions. 

7. Resolution to the Conjecture for f{z) having multiple singularities. 

For more general cases, we consider next the case of A'^g > 2 for singularities of two kinds, one 
being those of equal order and arbitrary locations whereas the other for those of arbitrary orders and 
locations, to begin with the case of A'^g = 2 for possible extension to multiple singularities. 
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7.1. Resolution for f{z) having two singularities of equal order. 

It is of interest to consider f{z) having multiple singularities of equal order, such as that found in 
the literature on water waves. In the case of A^s = 2 for two singularities of equal order, we have 

f{z) = Mi(zi-z)'= + M2(z2-z)'^ = S~o(^i^r" + i2 ^2"")7n(A;)(-zr, (19) 

where Li = MiZ^ {i = 1,2) and jnik) is given in (f7a|) . Similar to dS])-®, matching series (fT9ll for z = e*^ 
termwise with the prescribed series in ([5]) yields 

(Li zf " + L2 Z2-")7n = (-l)"cn (n = 0, 1, 2, • • • ), (20) 

These are the set of transcendental algebraic equations for solving the five base parameters in total, 
(Ml, M2, zi,Z2, k) and for verifying the consistency for this case. 

First, we eliminate the terms with L2 by operation for {{k — n)/{n + 1)^2 }cn + Cn+i, giving 

Li{z~^ - Z~^)'yn+lZi"' = (-l)"c„+i, Cn+l = , ^ , J! Cn + c„+i. (21a) 

{n+ 1)Z2 

Next, Li can now be eliminated by taking the ratio Cn+i/cn (n = 1, 2, • • • ), yielding 

= = —7— = —-rHnik, Z2) Hnik, z) = — — , (21b) 

Cn 7nZi n + lzi n + 1 nz - [n - 1 - k)/ Rn-i 

which can clearly be written in symmetry between the two singularities in equality as 

Gn{k,Zi) = Gn{k,Z2), Gn{k, z) = Hn{k, z) / z (n=l,2,---). (21c) 

Finally, expanding out (|21bp yields the basic equation for this case as 

{n + l)RnZiZ2-{n-k){{zi+Z2)-{n-k-l)/{nR^^i)} = {) (n=l,2,---), (22) 

which also exhibits the symmetry between the two singularities interchangeable in designation. 

To resolve the three unknown parameters {k, zi, Z2) in (j22p . we may take four leading equations of 
to eliminate {zi + Z2) and Z1Z2 in two steps, yielding one equation for A; (n = 1, 2, • • • ) as 

Ki-l,k) f Rn+i Rn+2 \ , K{l,k) i?„ , X(0,fc) Rn+2 _^ K{j k)^"^^^^ (23a) 



i?„_i \K{l,k) K{2,k)J Rn+i K{0,k) Rn K{2,k) ' ^ n + l+j' 

which is a cubic equation for k (after dismissing the complex conjugate roots by confirmation). With k so 
determined, (Z1Z2) can be deduced by eliminating zi + Z2 from giving for n = 1, 2, • • • , 

/ n — k n — /c — l\,/n+l„ n + 2 „ \ , , , ^ , , ^ 

\[n + l)Rn Rn-i J \n — k [n + 1) — k J 

for n — 1, say, with which zi + Z2 (= 2B{n,k), say) then follows from (|22p for n — I. Finally, zi and Z2 are 
obtained by combining ziZ2 — A(n, k) and zi + Z2 — 2B{n, k) to give the quadratic equation, 



2Bzi + A = Q — > zi = B+ ^B^ - A, Z2 = B - vW^^. (23c) 



At last, the solution is completed with the singularity magnitudes found from ()21ap for n = and verified 
for its feature of consistency using (jl9p for n > 5. Or in another approach, we may use (Mi, k, zi) Just 
obtained to provide its contribution, say cjj (n = 0, 1, • ■ ■ ) to differ from the original Fourier coefficients 
Cn by c„ — as the new data for the other singularity (M2, k, Z2) by direct application of ([8]) for verifying 
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the solution. Our resolution to the Conjecture is then accomplished for this case of f{z) having two 
singularities of equal order. When there are more than two singularities of equal order, this method can 
again provide a scheme for solution by induction. 

7.2. Resolution for f{z) having two individual singularities. 

For f{z) having two singularities of arbitrary orders and locations, we have 

/(z) = Mi(zi - zf' + M2(Z2 - zf^ = S^=o(ii^r"7n(A:i) + L2 z^''jr^{k2))i-zr, (24) 

where Li = Mizf (i = 1, 2). A gain, matching series (|24|) with series in (0) yields for this case 

^r"7n(A:i) + L2 ^2""7n(A:2) = (-Ifcn (n = 0,1,2,- ••), (25a) 

These are the set of equations for solving (Mi,Zi,ki, i = 1,2) and for verifications. Elimination of Li 
and L2 can be carried out in close analogy with that for (jl9p . yielding the basic equation as 

zi - Z2 - 1+ . (^2 - n) + (n + l)RnZ2 .25b) 

{k2 - n + l)zi - {ki-n + l)z2 ki-n + 1 nz2 - (n - /c2 - l)/-R„-i ' 

which reduces to (|22|) for ki = k2 = k. Pursuing the image symmetry that ()25bp is invariant if (zi, ki) 
and {z2, k2) are interchanged between the two independent singularities, we obtain 

Gn{ki,Zi) = Gn{k2,Z2), Gn{k, z) = Hn{k, z) / z (n = l,2,---), (26) 

where H{k^ z) is given in (|21b[) . 

Now, the four remaining unknowns, ki^k2,zi^ Z2 can be solved by taking four equations of (j25bp or 
(j26p . for n = 1,2,3,4, with the rest serving for verifying the solution. 

Example 2. Here we have the coefficients in numerics for series ([5]) as 

ci = 4, C2 = 1.5, C3 = 0.583333, C4 = 0.236111, C5 = 0.099537, cq = 0.043596, • ■ ■ , — > 
Ri = 0.375, R2 = 0.388889, R3 = 0.404762, R^ = 0.421569, R5 = 0.437985, • • • , 

Not knowing the number of singularities, we first test for A'^s = 1. Then, like in Example 1, we find 

k{l) = -1.074074, k{2) = -0.781818, A:(3) = -0.572727, • • • , 

for the singularity key index k{n) of the few leading orders. Their lack of any feature of consistency 
evidently indicates the singularity distribution being not in the single singularity group. 

So next, we test out the group of two separated singularities for this f{z). By substituting the above 
RnS in the four leading equations of (l25bp to solve for {ki, k2, zi, Z2), adopting successive elimination 
or any other adequate algorithms, we obtain, after some algebra, the final solution as 

ki = k2 = -1, zi = 2, Z2 = 3, Ml = -2, M2 = -9, 

in which Mi , M2 follow from ([25a]) for n = and n = 1 , with the consistency exhibited as 

112 

Gi{zi, k) = Gi{z2, k) = --, G2{zi,k) = G2{z2, k) = --, ^3(^1, k) = G3(z2, k) = --, 

followed by G^^zi, k) = G4{z2, k) = —5/6, each with a relative error of 10^^, and likewise for n > 5. 
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Example 3. As the two singularities in Example 2 turn out to be of equal order, we can use the same 
coefficients c„ for a comparative study. Substituting the RnS of Example 2 in (j23ap yields k = —1, 
with the other two complex conjugate roots dismissed by verification. With k = —1, ()23bp for n = 1 
produces ziZ2 = G and then ([22]) gives zi + Z2 = 5. Therefore by (|23cp . we obtain the solution 

k = -1, Zl = 2, Z2 = 3, 

in agreement with Example 2, here also with relative errors of 0(10^^). In comparison, the algebra 
involved here is noticeably simpler than the former case, especially for conducting verifications. 

When more than two singularities are encountered, this case for two singularities can be expected 
to provide a mathematical scheme for resolution by induction, i.e. by eliminating the magnitude of 
one singularity at a time till all the magnitude variables are eliminated for solving the locations of 
the singularities and their key order indices by complex algebra. However, the computation along this 
approach is still not simple, whereas a more general approach can now be addressed next. 

8. Resolution to the Conjecture involving multiple singularities by induction. 

For this general case, we reconsider the group of functions having multiple singularities, all lying 
outside the open domain P"*" bounded by the unit circle C, and situated all at different radial distances 
from the origin (as the case of singularities situated on a concentric circle beyond C has been discussed 
in §6). To begin with, we first consider f{z) having two singularities at zi and Z2, both real and positive 
(1 < Zl < Z2), a special case which is apt for illustrating the underlying basic principle more easily, and 
can be extended to functions having more singularities situated more freely. 

To have the singularities of f{z) determined by induction, we first consider the direct problem. 
8.1. The direct problem. The direct problem has the singularities all given in closed analytical form 
as a basis for developing a general method to determine the singularities by induction so that the method 
can readily be adapted to resolving the inverse problem as stated by the Conjecture. 

Thus, we take, parallel to (|24p . here with 1 < zi < Z2, 



f{z) = fi{z) + f2{z) = ^]=iMj{l - z/zjp = S~ Mj^n{kj){-z/zjr = SnC„z", 

cn = Mi7„(fei)(-zi)-" |l + m^^MA") (m = M2/M1, < A = zi/z2 < 1), (27) 



ln{k 

where Cn derives its value from the analytic formulas given for this direct problem to give 

^ ^ Cn+i ^ 1 - h/n On Q ^ 1 + m(7n+i(fc2)/7n+i(A:i))A"+^ ^^^^ 
Cn 1 + 1/n Zl' 1 + ?n(7„(/c2)/7n(A;i))A" 

For < A < 1, A" ^ and Qn ^ 1 as n ^ 00. Hence, for given A < 1, there exists N 3 Vn > 
N, |1 — Qn\ < e. Thus, as e ^ 0, Rn tends to that for a single singularity, which can then be 
determined in this asymptotic limit when the effects due to the singularity farther away fall off. 
Domb-Sykes plot. Proceeding for this asymptotic solution, it is apt to make the Domb-Sykes[ll] plot 
of Rn versus 1/n as a function i?„(l/n), just as so conspicuously exhibited by the functional relation in 
(j28p . with {ki, zi,Qn{X)) as the primary parameters. The basic principle and the main features of the 
solution can be illustrated below. 

Example 4. Function f{z) having two logarithmic singularities. Of basic interest is that 

f{z) = log{(l - Z/Zi)(l - Z/Z2f} = log(l - Z/Zi) + 21og(l - Z/Z2) (zi = 1.1, Zl < Z2), (29) 
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Here, we have f{z) = fi{z) + f2{z), ki = k2 = 0, m = 2, z\ = 1.1, A = z\l Z2 (0 < A < 1), hence 



n 



Qn 

n + 1 zi 



1 Qr. 



1 + 2A"+i 
1 + 2A" 



in 



1,2, 



(30) 



This function ii„(l/n) has been computed by employing Mathematica 7 for three values of A = 
0.125, 0.4, 0.8, over the range n = 1,2, ■•• ,100, with the resulting i?(l/n) plotted in three lines of 
sequence of points for each n as shown in Fig. 1, with the indicating the points for A = 0.125, the □ 
for A = 0.4, and the • for A = 0.8. Here, for n > A?^ ~ 20, the three lines of dotted points apparently 
merge into one straight line, ending at 1/n = 0.01 (or n = 100). This line is further extrapolated 
by being fitted to R = {po +pi/n)/{\ + qi/n) to reach the i?-axis intercept at i?(0) = 0.9091, of 
slope i?'(0) = —1 which gives A: = 0, thus implying fi{z) logarithmic and further zi being located at 
zi = l/i?(0) = 1.1, with almost no error. These are the two important data found from the plot. It is 



R(n) 
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Figure 1: Domb-Sykes plot of i?n(l/n) for three values of singularity radial-distance ratio A = ^1/2:2, shown by 
points for A = 0.125, by □ for A = 0.4, by • for A = 0.8; n = 1, • ■ • , 100. 



of interest to note that for n < N 20, the three dotted curves have the line (for Z2 = 8zi) staying 
the closest to the asymptotic line. La : Ra = R{0) + i?'(0)/n, whereas the □ line (for Z2 = 2.52;i) and 
the • line (for Z2 = 1.25zi) bifurcate increasingly more from the straight La line with decreasing radial 
distances. This signifies that the interaction between two point singularities increases with decreasing 
distances apart. This is another key feature demonstrated by the Domb-Sykes plot. 

At this point, we note that the number of singularities, Ng = 2, and k = being both known in this 
direct problem, we can directly solve for the other singularity by taking the first two Ri , R2 to deduce 
from their formulas to attain m = 2 and the values of Z2 for the three cases, and their magnitudes from 
Co to complete the solution to this direct problem. 

8.2. The inverse problem. Proceeding to resolve the Conjecture, we now take up the inverse 
problem which has only the Fourier coefficients Cn of ^ prescribed in numerics, ordinarily by the 
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solution to a problem solved by a numerical scheme, but here assumed to take the value c„ from the 
above direct problem for a single case for A = 0.4, say. With c„ = c„ known, we plot the coefficient 
ratio Rn = c„_|_i/c„ to attain in this case again the same result as that shown in Fig. 1, thereby having 
{ki = 0, zi = 1.1) determined for the first singularity of fi{z) (lying nearest to the unit circle C) as 

fiiz) = Ml log(l - z/zi) = c; - S„=i c^z" « = 0, 4 = Mizf "/n, n > 1), (31) 

where = Cq = 0, Mi = 1 by deduction from = Cn asymptotically, and hence fi{z) is all known. 

To this end, it is vital to note that the dotted line (here with □) being curved indicates that f{z) 
has more singularities, possibly not just one more, to complete for resolving the inverse problem. 
8.3. Resolving the inverse problem by induction. To continue, we rid f{z) of fi{z), giving 

f2{z)=f{z)-Mz) = ^oc:z^, C = cn-c!^ ^ < = (n = l,2,...). (32) 

With here all known, we can carry out the Domb-Sykes plot for R'^{l/n) to determine the primary 
parameters {k2, Z2) for f2{z)- In this particular case for f{z) of (j29p . we will find from the i?^(l/n) plot 
that the entire dotted line is straight, of slope dR" /d{l/n) = —1, and the /^''-intercept at l/z2 = 0.36364 
for Z2 = 2.75 = 2.5zi, as the last singularity of f{z), therefore completing this inverse problem. 

In general, for generic cases, the R"{\/n) line may appear again curved, in which case we repeat 
the analogous plot for R^'^\l/n)) for fs{z), • • • , until R^^\l/n)) for fe{z) bearing out a straight line in 
the plot for the very last singularity of fi{z) to complete the inverse problem. This then constitutes the 
general method of resolving the Conjecture by induction. 

9. Resolution to the Conjecture involving essential singularities. 

To enhance addressing this issue, we apply the generalized Hilbert transforms (cf. Wu[l]). 
9.1. Classical Hilbert transform. First, we consider a class of analytic function f{z) = 
u{x,y) + iv{x,y) which is regular in the upper half z-plane for Im z > 0, and vanishes as \z\ — > 00 
uniformly in < argz < vr, with the upper-half (or lower-half) z-plane designated as the open domain 
P"*" (or 'D~), bounded by the x-axis {y = 0) and the point z = 00. Then we have the Hilbert transform 
between u(x) and v{x) along the rr-axis given by 

n = H[v[x)] = -V / -, v{x) = H n ^ = — T> / , (33) 

with H denoting the Hilbert transform and the inverse transform. We note that as a point z tends 
to (x,0) from V+ (or V), f{z) in ([IaD-([Tb]) tends in the limit to /+(x) = /(x), /"(x) = 0, f{x) = 
u{x) + iv{x), u{x) and v{x) being given by (f33|) . all as a special case of ([3]) for this particular geometry. 
Example 5. Hilbert transform of sinusoidal functions. For v{x) = sinx, we have 

n(0 = H[sm{x)] = -V H ^J^^l^ = cos(e). (34) 

J -00 x-e, 

Hence, /(x) = u{x) + iv{x) = e*^ on the x-axis, which can be continued into the entire z-plane, giving 

/(x) = n(x) + iv{x) = e^^ — > f{z) = e''' (V z G < \z\ < 00). (35) 

Clearly, this function, e*^ = e*^~^ is regular in the upper-half z-plane, vanishing exponentially as 
y — s- +00, but becomes essentially singular as y ^ —00 in domain D~ (the lower-half z-plane). 
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9.2. Complement Hilbert transform. On the contrary, if complement function F{z) = U{x,y) + 

iV{x,y) is regular in the open domain 'D~{z\Im z < 0}, then applying the general formula ([3]) 
produces for this particular geometry in 'D~ the complement Hilbert transform between U{x) and V{x) 
along the x-axis the equations (cf. Wu[l]) as 

with H denoting the complement Hilbert transform and H ^ the inverse transform. They can be 
employed for evaluating transforms of various complement functions. 

Example 6. Complement Hilbert transform of sinusoidal functions. If for the complement 
function F{x) = U{x) + iV{x) on the x-axis, V{x) = sinx (equal in value to v{x) = sinx in (34), then 
by complement Hilbert transform (j36p . we obtain 

C/(0 = :^[sin(x)] =--V r ^^^^ = - cos(e). (37) 

vr y_oo x-i 

Hence by analytical continuation, we have 

F{x) = U{x) + iV{x) = -e"*^ — > F{z) = -e"*^ (V z G < |z| < oo). (38) 

Obviously, F{z) = — e"*^'"*"^ is regular in and has an essential singularity in D"^ at z = oo. Therefore, 
one equal value of one conjugate function, e.g. v{x) = V{x) = sinx, produces an analytic function f{z) 
regular in D"*" and another F(z) regular in 'D~ , whilst both having an essential singularity in each of 
their complement domain, respectively, in splendid analogy with the original Conjecture in §1 and the 
Complement Conjecture in §5, conjointly resulting in (I16p as their joint roots. 
10. Discussion and conclusion. 

In this study we have illustrated how the solutions can be obtained to the inverse problem for the 
generalized function f(z) having total number of A^^ = 1, 2, ■ ■ ■ singularities in resolving the Conjecture. 
The simpler cases (i)-(ii) for /(z) having only one or two singularities are seen to have provided the 
study a sound base by exhibiting a clear structure of all solutions to the inverse problem. The general 
method developed here in §8 for resolution to the Conjecture by induction is suitable for singularities 
located at differing radial distances so that the one closest to the origin can be singled out by applying 
the Domb-Sykes plot for an asymptotic solution, leaving the remaining singularities to be similarly 
resolved one at a time. This is complemented by the case of singularities lying on a concentric circle by 
the analytical method expounded in §7 and by another approach utilizing the series sign pattern criteria 
that has a longer history with the pioneering works of Van Dyke[15] and Li[14]. The success of these 
methods and the high accuracy of the general results are seen to stem from the central role played by 
the power series in complex form composed of orthogonal terms. These methods can now be employed 
for application and further development to all pertinent problems in mathematics, applied science, and 
mathematical physics. 

As mentioned in foregoing occasions, there are at least two major areas in analysis of holomorphic 
functions where the existing challenges can be benefited from the resolution of the Conjecture. The 
long and rich history of the fully nonlinear theory of dispersive water waves can serve as a splendid 
representative for all the disciplines in which the key advances may require a final resolution of their 
principal theory of perturbation expansions. In various other circumstances, new grasp of the exterior 
singularities in domain complement to the solution region of the primary variables may furnish crucial 
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data base to ascertain the certitude of solution validity and profound comprehension of the target phe- 
nomena. It may also cast light on the course of expounding the key mechanism underlying instabilities 
of a nonlinear system at an initial stationary state. 

In concluding, it is reasonable to say that this is just a beginning. New issues may arise to require 
answers. The exterior singularities of analytic functions may present a variety of new interests, e.g. 
coalescence of singularities with varying parameters, productions of essential singularities, etc. We may 
indeed expect to see more experience accumulating to enrich more of this new field. 
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